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We construct a TV— dimensional Gaussian landscape with multiscale, translation invariant, logarithmic 
correlations and investigate the statistical mechanics of a single particle in this environment. In the limit of 
high dimension N oo the free energy of the system in the thermodynamic limit coincides with the most 
general version of Derrida's Generalized Random Energy Model. The low-temperature behaviour depends 
essentially on the spectrum of length scales involved in the construction of the landscape. We argue that our 
construction is in fact valid in any finite spatial dimensions A'' > f . 
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The idea of energy landscapes pervades the theo- 
'retical description of glasses, disordered systems, pro- 
iteins, etc. [1]. The general goal is to classify typical 
'random potentials and establish their universal proper- 
ities, not unlike the Random Matrix Theory paradigm. 
This knowledge can then hopefully be used to describe 
'generic static and dynamic properties of complex sys- 
.terns, by addressing single point particle behavior in 
such potentials. In this respect, the Parisi solution for 
iSpin-glasses is fascinating: it reveals landscapes with 
'a surprisingly complex, hierarchical structure of valleys 
within valleys within valleys, etc. [5] . It is often however 
argued that the ultrametric properties of Parisi land- 
'Scapes are hardly compatible with a finite dimensional, 
^translation invariant space. 



particle, confined inside an A^— dimensional spherical 
box of radius L, is described by the coordinate vector 
r — (ri, ...,rAr), |r| < L. It feels a Gaussian-distributed 
random potential V{r) with zero mean, and with covari- 
ance chosen to be isotropic, translation invariant and 
with a well-defined large TV— limit: 
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In Eq.([T]) and henceforth the notation (. 
an ensemble average over the random potential, and / 
is a well behaved function of order unity. The thermo- 
dynamics of this model is described by the free energy: 
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In this paper we provide an explicit construction 
|of a Gaussian random potential in Euclidean, N di- 
'mensional spaces, with a specific form of long-ranged 
.correlations which reproduces all the features of Parisi 
'landscapes. More precisely, we show that the ther- 
modynamics of a single particle in a multiscale, loga- 
rithmically correlated potential is exactly described by 
Derrida's Generalized Random Energy Model (GREM, 
[3]), with an arbitrary (possibly infinite) number of lev- 
els of hierarchy. Although our proof concerns, strictly 
speaking, the limit TV — s- cx), we are confident that 
our results hold in arbitrary finite dimension > 1. 
This conviction is built both on physical arguments and 
on the beautiful results of Carpentier and Le Dous- 
sal [4] on the monoscale version of our model in fi- 
nite dimensions, which, as shown recently, match the 
exact results of the same model when ^ oo [5]. 
The model is defined as follows: the position of the 
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cxp —f3V{r) dr 
(2) 

as a function of the inverse temperature (3 = 1/T. 

Models of this kind has been studied extensively, and 
in the high-dimensional limit detailed analytical calcula- 
tions performed in [6| revealed that the nature of the low 
temperature phase is essentially dependent on the be- 
havior of the covariance f{u) at large distances. Namely, 
for short-ranged correlated potentials, the low tempera- 
ture phase turns out to be described by one-step replica 
symmetry breaking scheme of Parisi. In contrast, for 
the case of long-ranged correlated potentials with f{u) 
growing as u^^ , the full infinite- hierarchy replica sym- 
metry breaking (FRSB) scheme has to be used. 

The problem was reconsidered in much detail re- 
cently in 5J. In the limit N oo, one actually finds 
a true phase transition as a function of temperature 
provided the size of the confining sphere L is scaled 
as R\/N, with the parameter < i? < oo playing 
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the role of effective radius of the sample. As the ex- 
istence of a phase transition for finite-size systems is a 
kind of pathology of the infinite-dimensional approxi- 
mation, one is mainly interested in the thermodynamic 
limit R ~f oo. A simple analysis then reveals a special 
role played by logarithmic correlation function: 

f{u) = fo-9^\ii{u + a^), (3) 

where g, a and /o are given constants. In particular, 
only for such case the critical temperature Tc tends to 
a finite value: Tc{R oo) = g. Furthermore, the free- 
energy found in this limit is given by the well-known 
Random Energy Model expression J^. Interestingly, 
these results coincide precisely, up a trivial rescaling, 
with those obtained earlier for the same potential ^ at 
1 < < oo dimensions from a Renormalisation Group 
treatment In the standard interpretation, below Tc 
the partition function becomes dominated by a few sites 
with particularly low random potential, where the par- 
ticle ends up spending most of its time 8 . Note that 
the logarithmic growth of the variance of the potential 
appears naturally in various physical systems of diverse 
nature, see 01 H]. 

The main observation of the present paper is that 
the above picture, despite looking rather complete, still 
misses a rich class of possible behavior that survives in 
the thermodynamic limit i? — > oo. Namely, given any 
increasing positive function <&(?/) for < y < 1, we 
demonstrate below that if one considers potential cor- 
relation functions f{u) which take the following scaling 
form 

/(.) = -21ni?$(^l^i|i±p), 0<^.<i?^ (4) 

the thermodynamics of our system in the limit i? — > oo 
is precisely equivalent to that of celebrated Derrida's 
Generalized Random Energy Model (GREM)[3J. The 
REM- like case Eq.(I3]) turns out to be only a (rather 
marginal) representative of this class: ^{y) — g^y. 

Let us explain the motivation of the above form, 
which will make the physical interpretation of the re- 
sults quite transparent. The idea is to write V{y) as a 
(possibly infinite) sum of independent Gaussian poten- 
tials: V{y) = X^ili ^(i")) each with a covariance: 

{V, (ri) Vk {V2))v - S^.kNU (^(ri - v^f^ , (5) 

where fi{u) ~ —gf In (u + + af ) as in ([3]), but each 
with its own constant gi, and small-scale cutoffs cho- 
sen to grow as a power-law of the system size: at = R^^ 
with < z^i < 1. Taking the continuum limit K —t oo 



with a certain density p{v) of exponents Vi, we end up 
with: 

f{u) = - /" p{v)g^{v) In [u + a^ ^ R^") dv, < u < R^ . 
Jo 

(6) 

Now, introducing u + = R^^ and identifying with Eq. 
(|4|) in the R ^ oo limit, we obtain that the function <1> 
has the following representation: 

Hy)^y I p{i')9^{^)dv+ I i^pW)g\y)dv. (7) 

Jo Jy 

Note also that in this representation, $'(?/) = 
p{iy)g'^{v)diy > 0, and $"(?/) > 0, where the 
number of dashes here and below indicates the number 
of derivatives taken. The main result of this work is 
the following: depending on the nature of the spectrum 
of the exponents v, discrete or continuous, we will 
recover, in the thermodynamic limit, either the free 
energy of the original GREM with discrete hierarchical 
structure, or of its continuous hierarchy analogue (see 
(ITS)) below) analysed recently in much detail by Bovier 
and Kurkova [TU], see also [TT]. 

The physical interpretation of our results is as fol- 
lows. Instead of one localisation transition tempera- 
ture Tc at which the particle chooses a finite number 
of "blobs" of size 0{a) where the potential is particu- 
larly deep, there appears K different transition temper- 
atures, where the particle localizes on finer and finer 
length-scales. The largest transition temperature Ti 
corresponds a condensation of the Boltzmann weight 
inside a few blobs of large size 0{R), but the particle 
is still completely delocalized inside each blob. As the 
temperature is reduced, the REM condensation takes 
place over smaller blobs of size 0{R^) inside each al- 
ready occupied large blobs, and this scenario repeats 
itself as the temperature is reduced, each time "zoom- 
ing" in on a smaller scale [T^] . 

The equilibrium free energy per degree of freedom 
of our model, Foo = limjv^oo -Fjv/A^ with defined 
in Eq. ([2]), can be found in a standard way for any 
covariance f{u) using the replica trick. The details of 
the corresponding analysis can be found in 0, and we 
give below a summary of the most essential formulae 
for the FRSB situation. For finite i?, the low temper- 
ature phase is characterised by the existence of a non- 
trivial, non-decreasing function cc(g), q S [go,9fc], with 
the two parameters go and qk satisfying the inequality 
1^ qa ^ qk ^ qd ^ ■ The corresponding Foo can 
be written in terms of only those two parameters, see 
Eq.(58) of [SI. Here we choose instead to introduce, 
along the line of the physical discussion given above, 
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two characteristic "blob" sizes (actually size squared) 
rfmin = - qk, rfmax = - 9o in tcrms of which: 



T. 



--ln[2.ed_] , ^ 



/(O) - dmin/'(dmin)] 



(8) 
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where dmin ^ c?max can be found for a given temperature 
T from the equations 



< dmin = 
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(9) 

Finally, the Parisi order-parameter function, which 
takes the values between and 1 and is the main mea- 
sure of the ultrametricity in the phase space, has the 
following shape 



xid) 



T f"'{d) 

' 2 [/"(d)]3/2 



yd e [rfmin, rfmax] 



(10) 



where we performed the overall change q ^ d = B? — q 
in comparison with [5] . This function must be now non- 
increasing, and one can verify that this is precisely the 
case e.g. for the family f{u) in Eq.(j4]). 

The above solution is valid for the temperature range 
< r < Tc, where the critical temperature is given 
in terms of the largest blob size dmax as: 



Fq ^max yf (^max) ■ 



(11) 



Above this temperature the solution is replica- 
symmetric (RS), corresponding to a delocalized phase 
for the particle: no particular region dominates the 
partition function. The corresponding free energy is 
given by: 



p — 



T 



In [27r(is 
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where satisfies 



ds=R' 



y2 



(12) 



(13) 



We now consider specifically correlation functions 
/(u) of the form ([4]). In what follows we will use the con- 
venient notations z = (21ni?)~^ and y — zln(ii + a^). 
As noted above, our multiscale logarithmic model en- 
sures that $'(?/) > and $"(j/) > for any < y < 1. 
We will assume for simplicity $'(0) — 0, relegating con- 
sideration of the general case to an extended publication 
[TB] . We start our analysis assuming the function <&"(y) 
is finite and differentiable, but later on will relax this 
condition. Our first goal is to find the largest blob size 



c^max from second equation in Eq.Q, and then to de- 
termine the critical temperature Tc. Introducing the 
scaling variable 2/max = ^In (dmax + a^), in the thermo- 
dynamic limit z — > we can look for a solution Uraaxiz) 
as a power series of z. One immediately checks that 
Vmayi{z) = 1 — z In 2 + O(z^). This imphes that the 
largest blob size is of the order of the system radius: 
rfmax ~ i?^/2 > for R ^ oo. F,q.^ then yields 
the critical temperature given in the thermodynamic 
limit by a very simple expression Tc — y'^P^l). Phys- 
ically, at Tc, the sample breaks up into blobs of size 
o(i?) and only a finite number of these blobs are visited 
by the particle. However, within each blob, all sites 
are more or less equivalent. Now we can treat along 
the same lines the first equation in ^ to determine 
the smallest blob size dmin for T < Tc- It can again 
be conveniently written in terms of the scaling variable 
ymin Z ln(dmin + a^)- In the thermodynamic limit 
z — > 0, it is again natural to look for a solution ?/min as 
a power series of z, in which we only retain the first two 
terms: ?/niin = iy^,+cz+0{z'^). Due to our assumption on 
differentiability of the function $' (y) we expand around 
y = ly,:, and after a simple calculation find c — I. This 
means that dmin behaves like c?min = ei?^"* for R oo, 
where i^, satisfies the equation 



r2 ^ 



(14) 



Since the function ^'{y) is monotonously increasing for 
y > 0, and $'(0) = we find that in the limit R ^ oo 
(i.e. z ^ 0), the equation Eq. lH)) must have a unique 
solution < v*{T) < 1 in the range of temperatures 
< T < = V*'(l)- In this regime, dmin < dmax- 
Physically, sites within blobs of size dmin or smaller are 
not resolved by the particle, which visits all of them 
more or less equally. 

Now we can easily find the free energy Fco by substi- 
tuting these results to Eq. ((8]) and extracting the leading 
term in the thermodynamic limit z ^ 0. We find the 
equilibrium free energy to be of the form Foo — (In R) T , 
where for < T < 

-^^TMT) + M^^ + 2/;Vl^d,. 

(15) 

For T > Tc the solution of ([13]) in the hmit i? ^ oo is 
given by dg = R^ t^+t'^ and substituting this to we 
find that the free energy is given by: 



-T = T 



[$(1) - $(0)] 



(16) 



Last but not least, we can determine the thermody- 
namic limit of the order-parameter function x{d) given 
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by Eq. (jlOp . which determines in a precise way how the 
particle locahzes on different scales. To leading or- 
der in z we find /"(u) = $'(y)/(w + a^f, f"'{u) = 
— 2$'(2/)/(m + a^)'^ with y — zln(w + a^). Introduc- 
ing again the scaling variable u = ^ ^ 
[ei?^''* , i?^/2] we see that the order-parameter function 
assumes the limiting form: 

This completes our solution of the problem for the 
case of continuous function ^'{y). At this point it is 
rather informative to consider the case of a discrete 
spectrum of K exponents i^i, i = satisfying 
^ < vk < vk-\ < ■ ■ ■ < v\ < ^ \. This corresponds 
to K superimposed logarithmic potentials with 

K 

g\v)p{v)=Y.3l^{^-^^\ (18) 

1=1 

with 5(u) standing for the Dirac delta-functions. 
The corresponding $'(y) consists of steps: 
^'(2/) — 9l^{v~^i)- A simple consideration shows 

that our earlier analysis for the values of dmax and the 
critical temperature still hold for such a case, so 
d„,ax = i?V2, and = [$'(!)] ^ ^ ai+ gi^ . . .g\,. 
The equation ^ used to determine dmin = i?^''™'- — 
now takes the following form: 

, , 1 - a^e-?'™/^ 1 , , 

T^=>gf 7 7^, z = . (19 

i—l 

A little thought shows that the solution should al- 
ways be in the form ymin — Vp -\- CpZ for small z, 
where the index p runs successively through the val- 
ues 1,...,K when decreasing temperature from Tc to- 
wards T — 0. Introducing a decreasing sequence of 

characteristic temperatures Tp — \jY^f^p gf, we find 
in the interval Tp+i < T < Tp the value ymin = 
i/p-Fzln - T^+i)/{T^ - T2). Thus, the value of 2/,„in 
jumps (and thus the size of the smallest frozen blobs 
c^min) when crossing each of the temperatures Tp, with 
the highest one being Ti — Tc. It is also clear that for 
Up < V < Vp-i one finds: x(v) = T/Tp when T < Tp. 
Since Tp and Vp decrease as p increases, x{v) for a given 
temperature T < Tc is step-wise constant with jumps 
at each Vp; the smaller u (i.e. the smaller the size of the 
blobs), the larger x{i'), meaning that the localisation 
effect is weaker and finally disappears when x{i') —>■ 1. 



The expressions for (2/min, 2/max) suffice to calculate 
the free energy expression in the thermodynamic limit. 
In the temperature range Tp+i < T < Tp we have 

P ^ K 

-T^Tvp + 2Y,{v,-i-Vi)T + - J2 [i^r-i- Vi)Tf, 

i—l 2— p+1 

(20) 

whereas for T > Ti = Tc the RS expression is 
1 

-^ = T+-Y,i^.,,^-u.)T^. (21) 

Interestingly, these expressions reproduce exactly, mu- 
tatis mutandis the free-energy of Derrida's GREM [3l 
110] ■ with a particularly clear interpretation in terms of 
particle localization inside smaller and smaller blobs as 
the temperature is reduced. 

Remembering the mentioned perfect match between 
the results of [3] and [5] in the limiting case Eq.([3]) 
it is very tempting to conjecture that the GREM be- 
haviour revealed by us in the infinite-dimensional set- 
ting should also hold in all spatial dimensions, down to 
N = 1, albeit with the largest exponent < 1. In- 
deed, essentially the same mechanisms are at play in 
both situations. We hope that the corresponding RG 
and travelling wave formalism of [4] can be generalized 
to support this conclusion. For finite values of K, where 
lengthscales are well separated, this looks indeed quite 
feasible. 

If this conjecture is true, we would then have in- 
deed explicitly constructed a Parisi landscape in finite 
dimensions fully in terms of stationary Gaussian pro- 
cesses. How do we reconcile this with the ultrametric 
properties of the Parisi construction? Consider the fol- 
lowing distance Dji defined for any two points r, r' in- 
side a sphere of the radius R in the Euclidean space of 
any dimension: 

Inflr-r'P + a^l 
^^^'■''■^ ^ 21n!R ' ^ < ' - ^ ^^^^ 

Parameterizing |r| = R°'^^\ < a < 1, we see that 
in fact limj^^oo ^'flXr, r') = max{Q!(r), Q;(r')} . The lat- 
ter function used as a distance converts the Euclidean 
sphere into a so-called ultrametric space: every triangle 
will have at least two sides equal. We thus conclude 
that in our model the covariance of the random poten- 
tial depends only on the ultrametric distance inside our 
growing sphere, not unlike the original construction of 
GREM [31 [TU] or directed polymers on a tree with disor- 
dered potential [14] , cf. discussions in [9l [4] in the single 
scale case. 
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Several aspects of the model deserve in our opin- 
ion further investigations, some of them to be discussed 
elsewhere |13) . In particular, the rich behaviour found 
in the thermodynamics of a single particle should also 
have interesting dynamical counterparts, cf. [iSl [16]. 
One also can study multifractality exponents reflecting 
the spatial organization of the Gibbs-Boltzmann weights 
and the associated singularity spectrum. Finally, let 
us mention that in = 1 the monoscale logarithmic 
landscape model has in fact deep connections with the 
muZtifractal Random Walk construction suggested in 
[17j . The present model suggests a natural generali- 
sation to a multiscale logarithmic processes [13] . 
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